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A feasibility study is presented for the search of the lightest bottom squark (sbottom) in a com¬ 
pressed scenario, where its mass difference from the lightest neutralino is 5 GeV. Two separate 
studies are performed: (1) final state containing two VBF-like tagging jets, missing transverse en¬ 
ergy, and zero or one 5-tagged jet; and (2) final state consisting of initial state radiation (ISR) jet, 
missing transverse energy, and at least one 5-tagged jet. An analysis of the shape of the missing 
transverse energy distribution for signal and background is performed in each case, leading to signif¬ 
icant improvement over a cut and count analysis, especially after incorporating the consideration of 
systematics and pileup. The shape analysis in the VBF-like tagging jet study leads to a 3cr exclusion 
potential of sbottoms with mass up to 530 (462) GeV for an integrated luminosity of 300 fb -1 at 14 
TeV, with 5% systematics and PU = 0 (50). 

Introduction - Weak-scale supersymmetry addresses 
the hierarchy problem, gives gauge coupling unification, 
and (in imparity conserving models) provides a robust 
dark matter (DM) candidate, the lightest neutralino 
(x5)- As such, it is one of the most widely studied frame¬ 
works for physics beyond the Standard Model (SM). 

The exclusion bounds on supersymmetric colored par¬ 
ticles belonging to the first two generations are already 
quite strong. For comparable squark (q) and gluino (g) 
masses, the data eliminates these particles up to approx¬ 
imately 1.5 TeV at 95% C.L. with 20 fb -1 of integrated 
luminosity m ■ 

On the other hand, the bounds on the masses of the 
colored third generation are much weaker due to smaller 
production cross section. Given that a new boson con¬ 
sistent with the SM-like Higgs has been observed, with 
mass in the region of 125 GeV [6], a weakly coupled light 
scalar must have its mass stabilized against quantum cor¬ 
rections. Probing top squarks (stops) is a high-priority 
study for the future given its importance in this context. 

Since left handed bottom squark (sbottoms, b) and 
stops come in the same electroweak doublet, it is equally 
important to search for light sbottoms. Moreover, light 
sbottoms can play a role in obtaining the correct relic 
density of a neutralino DM candidate, through coannihi¬ 
lation effects Sbottom pairs produced from QCD in¬ 
teractions will decay to the lightest superpartner, which 
we will assume to be a stable neutralino, through the 
process b fry?. When the mass difference between the 
sbottom and the neutralino is large, the 6-tagged jet is 
sufficiently boosted. In that case, the standard procedure 
for the sbottom search involves final states containing two 
jets, at least one of which is 6-tagged, and large missing 
transverse energy p? coming from the neutralino. This 


has been the strategy of sbottom searches at CMS and 
ATLAS. 

The current exclusion bounds on sbottoms are as fol¬ 
lows. With 20 fb -1 of data at 8 TeV, the ATLAS Collab¬ 
oration has ruled out sbottoms up to 650 GeV, for neu- 
tralino masses less than 300 GeV [8], and up to 255 GeV 
for the mass-degenerate scenario [22]. Similar exclusion 
bounds have been obtained by the CMS Collaboration 
when the similar search strategy is employed mm ■ 

The purpose of this paper is to propose a new search 
strategy for sbottoms in compressed regions of parameter 
space, where the m^ — m^o mass difference is small (we 
have kept this value at 5 GeV throughout this study). 
The challenge in the compressed region is that both p? 
and the px of the 6-tagged jets become small, due to in¬ 
sufficient boosting of the objects coming from the sbot¬ 
tom decay. We point out that this can be overcome in 
final state topologies containing two boosted forward jets 
in opposite hemispheres (reminiscent of vector boson fu¬ 
sion jets). Gluons radiated off of the forward jets can 
produce sbottoms which decay to 6-tagged jets and pr 
in the central region of the detector. To balance the high 
initial px of the incoming partons, the centrally produced 
decay products are boosted, even in the compressed re¬ 
gion. 

This topology has been proposed by some of the au¬ 
thors as a probe of both colored and non-colored super- 
symmetric particles. Charged and neutral Wino produc¬ 
tion followed by decays to y^ y i a a light slepton has been 
studied in muni, while searches for Winos and Higgsi- 
nos in the final state of two VBF-tagged jets and pr has 
been proposed in m- Moreover, top squarks have been 
studied in the compressed region using this topology [12]. 

The search strategy using the initial state radiation 
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(ISR) jet has been already employed by the CMS and 
ATLAS Collaborations, and has shown a good sensitiv¬ 
ity for signals with compressed spectra tang. We also 
consider a separate study with final state consisting of 
an ISR jet, pr, and at least one 6-tagged jet. 

In each case, the shapes of the pr distributions for sig¬ 
nal and background are studied using a binned likelihood 
following the test statistic based on the profile likelihood 
ratio. We find that shape analysis, compared to a simple 
cut and count analysis, yields a significant improvement 
of the compressed sbottom mass reach after incorporat¬ 
ing the effects of systematics and pileup. 

In the rest of the paper, we present results first from 
the analysis with VBF-like tagged jets in the final state, 
and next the analysis with a ISR jet in the final state. 
We end with our discussions. 

VBF-like Tagging Jets Study - For this feasibil¬ 
ity study, inclusive bb* + multijets samples are gener¬ 
ated with b masses in the range of 15-1000 GeV, keeping 
m b ~ m xi ~ ^ GeV. Both QCD and weak production 
processes are included. The Xi in our studies is mostly 
Bino and the sbottom decays entirely through the canon¬ 
ical channel b —>> bx\. The other colored particles, neu- 
tralinos and charginos are assumed to be much heavier. 

Signal and background samples are generated with 
MADGRAPH5 [13] followed by the parton showering and 
hadronization with PYTHIA M and the detector simula¬ 
tion using DELPHES m One advantage of the DELPHES 
simulation is that it can simulate pileup pp interactions, 
which was not possible with other fast simulation pro¬ 
grams available in the HEP community. We used the 
Snowmass detector configuration as defined in m , which 
represents the typical performance of the CMS and AT¬ 
LAS detectors. We perform a 14 TeV study for pileup 
PU = 0, 50, and 140 scenarios, with an assumption of 
5% systematics on the signal and background. 

Given the compressed spectrum, we explored the sen¬ 
sitivity with two VBF-like tagged jets and large missing 
transverse energy final state. We also found that with the 
VBF-like event signature, the b jets from the sbottom de¬ 
cays also get boosted in the central part of the detector, 
providing another approach for probing this compressed 
signal. Overall, two studies are performed: (1) final state 
with two VBF-like tagged jets, pr and zero 6-tagged jet 
and (2) final state with two VBF-like tagged jets, pr, 
and one b -tagged jet. 

The cuts employed are as follows. 

(1) pr~pr asymmetry cut: the condition 


\Pt ~ Pr\ 

pT + pT 


< 0.2 (0.5) for PU = 0(140) 


(1) 


is imposed to protect against occasional loss of high pr 
jets due to the aggressive pileup subtraction in DELPHES . 
Here, pr i s defined as the negative vectorial sum of jets 
with px > 30 GeV, muons, electrons, and photons. The 
p T was found to less pileup dependant, and used instead 
of pr in this analysis. 


(2) Boosted jet cuts: the event is required to have 
the presence of at least two jets (ji, j 2 ) satisfying: (i) 
pr(j l) > 50(200) GeV and p T (j 2 ) > 50(100) GeV for 
PU = 0(140) in \rj\ < 5; (ii) \Arj(j 1 J 2 )\ > 4.2; (Hi) 
rijirjj 2 < 0; (iv) dijet invariant mass Mj 1 j 2 > 1500 GeV; 
(v) missing transverse energy p T > 50 GeV. 

(3) Vetoes: We veto electrons, muons and tau-tagged 
jets. In the final state study with two VBF-like tagged 
jets, pr and zero 6-tagged jet, a 6-tagged jet veto is also 
applied at this stage. 

(4) 6-tagged jet requirements (for 1 6-tagged jet study): 
For the tagged jets plus pr plus b study, we require 
exactly one tight 6-tagged jet in the final state, with 
Px < 80 GeV to suppress the tt background. 

(5) p T cut: We require pr > 200 GeV. 

(6) A <pjj\ From the search for invisible Higgs boson in 
the vector boson fusion at CMS [20], the QCD back¬ 
ground is reduced to a low level by requiring the az¬ 
imuthal separation between the VBF-tagged jets to be 
small. Here, we require A (f)jj < 1.8. 

The cut flow table for the benchmark point for PU 
= 0 in the VBF-like tagged jets, pr, and zero 6-tagged 
jet study is presented in Table [I| 


TABLE I: Cut flow table for m g = 500 GeV with PU = 0, in 
the final state with two VBF-like tagged jets, pr, and zero 
6-tagged jets. 


Selection 

Signal (pb) 

Background (pb) 

Boosted jets 

5.6 -10“ 3 

10.2 

6-tagged jet veto 

5.5 -10“ 3 

9.8 

Lepton veto 

5.4 -10“ 3 

6.5 

pr > 200 

3.3 10“ 3 

0.6 

A 0 

2.1 -10 -3 

0.25 


Figure [l] shows the distributions of pr normalized to 
unity for signal (green dotted histogram) and the dom¬ 
inant V(W, Z)+jets background (black solid histogram) 
after all selections except pr requirement, for the bench¬ 
mark point with m^ = 500 GeV, m^o = 495 GeV, in 
the case of PU = 0 for the VBF-like tagged jets plus 
pr study. Based on this distribution, and similar p T 
distributions for the VBF-like tagged jets, pr , plus one 
6-tagged jet study, a shape analysis was performed with 
different pileup scenarios. A local p-value is calculated 
as the probability under a background only hypothesis 
to obtain a value of the test statistic as large as that 
obtained with a signal plus background hypothesis. The 
significance z is then determined as the value at which 
the integral of a Gaussian between z and oo results in a 
value equal to the local p-value. 

In Figure [2j we show the significances of the com¬ 
pressed scenario with — m^o = 5 GeV, at 300 fb -1 
with the cut and count method (left panel) and the shape 
analysis (right panel), using joint likelihood to combine 









3 



FIG. 1: Distributions of If? normalized to unity for sig¬ 
nal (green dotted histogram) and dominant V+jets back¬ 
ground (black solid histogram) after all selections except Jjh 
requirement, for the benchmark point with = 500 GeV, 
rri-o =? 495 GeV, in the case of PU = 0, for the channel with 
VBF-like tagged jets, J]L r, and zero 6-tagged jet. 

the studies with and without 6-tagged jets. A systematic 
uncertainty of 5% is uniformly assumed. From top to bot¬ 
tom, the black solid, red dashed, and green dotted curves 
show the cases of PU = 0, 50, and 140, respectively. The 
red solid horizontal lines denote the 3cr and 1.69a levels. 
The shape analysis leads to a 3a exclusion potential of 
sbottoms with mass up to 530 (462) GeV with 5% sys- 
tematics and PU = 0 (50). (The 95% CL exclusion reach 
for the PU = 50 case at 300 fb -1 is 541 GeV.) For the 
most conservative PU = 140 case, using the shape (cut 
and count) analysis, it is possible to probe compressed 
sbottoms at the 3a level up to m~ h — 380 (300) GeV. 

ISR Monojet study - We now turn to our second 
analysis, in which the final state consists of an initial 
state radiation (ISR) jet, If t, and at least one 6-tagged 
jet. Our analysis mostly follows m , which studied this 
scenario with selections optimized for the 8 TeV LHC. 
The event selection is as follows. The If x asymmetry cut 
is applied as above to protect against occasional loss of 
high pt jets due to pileup subtraction. The leading jet is 
required to be non 6-tagged, and have px > 120 GeV. At 
least one 6-tagged jet with Pt > 25 GeV and |^| < 2.5 
is required. The leading 6-tagged jet is required to satisfy 
Px(5i) < 100 GeV, with A0(p T (5i), ^x) < 1.8. Leptons 
are vetoed. We also require If x > 430 GeV, and the top 
quark transverse mass to satisfy Mf > 200 GeV. 

The If T distribution (after all cuts except If' x cut), nor¬ 
malized to unity, for signal (red dotted histogram) and 
dominant V+jets background (black solid histogram) is 
shown in Figure [3l for m^ = 500 GeV and m^o = 495 
GeV, in the case of PU = 0. 

A shape analysis is performed following the method 
described above. This yields the significance plot shown 
in Figure 0 for m= 5 GeV with 300 fb 1 of 
data. From top to bottom, the black solid, red dashed, 
and green dotted lines show the cases of PU = 0, 50, and 


140, respectively. The red solid horizontal lines denote 
the 3a and 1.69a levels. In the ISR + 6-tagged jets study, 
we found the 3a level reach to be 250 GeV for the PU 
= 0 case. 

Discussion of Snowmass Simulation- For the 

study with VBF-like tagged jets, we have performed a 
joint analysis of the zero and one 6-tagged jet final states 
and displayed the exclusion reach in the right panel of 
Figure [2j Figure [5] shows the separate significances in 
the two final states, after shape analysis, assuming PU 
= 0 and a systematic uncertainty of 5%. It is clear that 
the performance is dominated by the zero 6-tagged jet 
channel, due to the small mass difference of 5 GeV be¬ 
tween the sbottom and Xi • For larger mass difference, it 
is expected that the final state containing one 6-tagged 
jet would perform better. 

In the one 6-tagged jet analysis, we found that the 
6-tagging efficiency at low px is critical in probing com¬ 
pressed scenarios. From the Snowmass detector simu¬ 
lation, the efficiency of 6-tagging reach 60% at jet px 
around 100 GeV [16], while the CMS and ATLAS have 
showed the ability to tag 6 jet with 60% efficiency for 
jets with p T above 30 GeV [17] 23]. Thus the VBF-like 
tagged jet plus one 6-tagged jet analysis can be signifi¬ 
cantly improved with more efficient and robust 6-tagging. 

Figures [2] and [4] shows the degradation in significance 
with higher number of pileup for both VBF-like and 
Monojet analysis. From [24], the expected jet perfor¬ 
mance of the CMS detector in pileup of 50 is compa¬ 
rable with the performance simulated in the Snowmass 
sample [16], while the Snowmass samples have much de¬ 
graded performance in the 140 pileup condition. With 
the upgraded CMS and ATLAS detectors optimized for 
the high-luminosity LHC condition and development in 
the pileup mitigation technique, we are expecting bet¬ 
ter physics object performance in the 140 pileup scenario 
and thus improved reach of sbottoms in the 140 pileup 
condition at HL-LHC. 

Conclusions - The main result of this paper is that 
the boosted jet topology can provide a feasible strategy 
to search for compressed bottom squarks with mass dif¬ 
ference of sbottom and lightest neutralino being 5 GeV. A 
shape based analysis is used to estimate the significances. 
There is 3a exclusion potential up to 530 (462) GeV for 
an integrated luminosity of 300 fb -1 at 14 TeV, with 5% 
systematics and PU = 0(50). We also performed an ISR 
+ 6-tagged jet study, and found the exclusion reach to 
be 250 GeV for PU = 0. 
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FIG. 2: The exclusion reach for the compressed sbottom scenario with mi — o = 5 GeV. A systematic uncertainty of 5% 

is assumed throughout. Left panel: The significance at 300 fb -1 as a function of m g in the cut and count method, for the 
joint studies with and without 5-tagged jets. From top to bottom, the black solid, red dashed, and green dotted lines show the 
cases of PU = 0, 50, and 140, respectively. The red solid horizontal lines denote the 3 a and 1.69cr levels. Right panel: The 
significance at 300 fb -1 as a function of mi with the shape analysis method, for the joint studies with and without 5-tagged 
jets. The legend is identical to the left panel. 



FIG. 3: ISR + 6-tagged jet study: Distributions of J]L t 
normalized to unity for signal (red dotted histogram) and 
dominant V+jets background (black solid histogram) after 
all selections except ljL r requirement, for the benchmark point 
with mi — 500 GeV, ra~o = 495 GeV, in the case of PU = 0, 
for the channel with ISR jet, one 5-tagged jet, and mht. 
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FIG. 4: ISR + 5-tagged jet study: The exclusion reach is 
shown for the compressed sbottom scenario with m g — m ^.o = 
5 GeV, for the channel with ISR jet, at least one 5-tagged jet, 
and mht. The significance at 300 fb -1 as a function of m g 
after shape analysis is displayed. A systematic uncertainty 
of 5% is assumed. From top to bottom, the black solid, red 
dashed, and green dotted lines show the cases of PU = 0, 50, 
and 140, respectively. The red solid horizontal lines denote 
the 3<7 and 1.69a levels. 
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FIG. 5: VBF-like tagging jets study: A comparison of 
the exclusion reach with m g — m^o m 5 GeV, for final states 
containing zero 5-tagged jet (red dotted curve) and one b- 
tagged jet (black solid line), along with the VBF-like tagged 
jets and The significance at 300 fb -1 is given as a function 
of mi after shape analysis, assuming PU = 0 and a systematic 
uncertainty of 5%. 










